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ABSTRACT 

Cylindrical algebraic decomposition is one of the most im- 
portant tools for computing with semi-algebraic sets, while 
triangular decomposition is among the most important ap- 
proaches for manipulating constructible sets. In this paper, 
for an arbitrary finite set F C R[j/i, . . . ,y n ] we apply com- 
prehensive triangular decomposition in order to obtain an 
^-invariant cylindrical decomposition of the n-dimensional 
complex space, from which we extract an F-invariant cylin- 
drical algebraic decomposition of the n-dimensional real spa- 
ce. We report on an implementation of this new approach 
for constructing cylindrical algebraic decompositions. 

1. INTRODUCTION 

Cylindrical algebraic decomposition (cad) is a fundamental 
and powerful tool in real algebraic geometry. The original al- 
gorithm introduced by Collins in 1973 [TTJ has been followed 
by many substantial ameliorations, including adjacency and 
clustering techniques [2], improved projection methods [25] 
CEHirjSE], partially built CADs [13] [26l [29], improved stack 
construction [14] and efficient projection orders [16] . 

The main application of CAD is quantifier elimination (QE) 
for which other approaches are also available. Some of them 
have more attractive complexity results [4] than CAD. How- 
ever, as pointed out by Brown and Davenport in [8], "there 
is the issue of whether the asymptotic cross-over points be- 
tween CAD and those other QE algorithms actually occur in 
the range of problems that are even close to accessible with 
current machines". In addition, these authors observe that 
CAD can help solving certain QE problems [7] [T5] that other 
QE algorithms can not. 



For a finite set F n C R[i/i, . . . , y n ] the CAD algorithm [TT] 
decomposes the real n-dimensional space into disjoint cells 
Ci, . . . , C e together with one sample point Si € Ci, for all 
1 < i < e, such that the sign of each f £ F n does not change 
in d and can be determined at Si. Besides, this decompo- 
sition is cylindrical in the following sense: For all 1 < j < n 
the projections on the first j coordinates (yi, . . . , t/j) of any 
two cells are either disjoint or equal. We will make use of 
this notion of "cylindrical" decomposition in C n . 

The algorithm of Collins is based on a projection and lift- 
ing procedure which computes from F n a finite set F n -i C 
R[yi, . . . , 2/n-i] such that an J^-invariant CAD of R" can be 
constructed from an f^-i-invariant CAD of R n_1 . This con- 
struction and the base case n = 1 rely on real root isolation 
of univariate polynomials. 

In this paper, we propose a different approach for comput- 
ing CAD, which proceeds by successive transformation of an 
initial decomposition of the complex n-dimensional space. 
Our algorithm consists of three main steps: 

Initial Partition: we decompose C n into disjoint constructible 
sets Ci, . . . , C e such that for all 1 < i < e, for each 
/ 6 F n either / is identically zero in d or / vanishes 
at no points of d. 

Make Cylindrical: we transform the initial partition and ob- 
tain another decomposition of C" into disjoint con- 
structible sets such that this second decomposition is 
cylindrical in the above sense. 

Make Semi-Algebraic: from the previous decomposition we 
produce an _F n -invariant CAD of R". 



Our first motivation is to understand the relation and possi- 
ble interaction between cylindrical algebraic decompositions 
and triangular decompositions of polynomial systems. This 
latter kind of decompositions have been intensively studied 
since the work of Wu [32]. The papers [3] [20] and book [30] 
contain surveys of the subject. The primary goal of trian- 
gular decompositions is to provide unmixed decompositions 



1 



of algebraic varieties. However, the third and fourth au- 
thors have initiated the use of triangular decompositions in 
real algebraic geometry |35| . Moreover, real root isolation 
of zero-dimensional polynomial systems can be achieved via 
triangular decompositions [331 1341 110] . 

A second motivation of this work is to investigate the possi- 
bility of improving the practical efficiency of CAD implemen- 
tation by means of modular methods and fast polynomial 
arithmetic. Such techniques have been successfully intro- 
duced into triangular decomposition methods [151 1241 122] . 
Each of the three main steps of the algorithm proposed in 
this paper relies on existing sub- algorithms for triangular 
decompositions taken from 28, 9, 34 b and for which efficient 
implementation in the Regular/Chains library [5T] is work in 
progress based on the highly optimized low-level routines of 
the Modpn library [2"3"] . 

Our third motivation is to extend to real algebraic geome- 
try the concept of Comprehensive Triangular Decomposition 
(ctd) introduced in [9]. The relation between CAD and 
parametric polynomial system solving is natural as pointed 
in |17] and the presentation therein of Weispfenning's ap- 
proach [3T] for QE based on comprehensive Grobner bases. 
This suggests that the algorithm proposed in this paper 
could support a similar QE method. 

This paper is organized as follows. A summary of the theory 
of triangular decomposition is given in Section [5] Section [3] 
and Section [4] are dedicated to the first two main steps of 
our algorithm whereas Sections [5] presents the last one. In 
Section[S]we report on a preliminary experimentation of our 
new algorithm. No modular methods or fast polynomial 
arithmetic are being used yet and our code is just high- 
level Maple interpreted code. However our code can al- 
ready process well-known examples from the literature. We 
also analyze the performances of the different main steps 
and subroutines of our algorithm and implementation. This 
suggests that there is a large potential for improvement by 
means of modular methods, for instance for the computation 
of GCDs, resultants (and the discriminants) of polynomials 
modulo regular chains. 

2. TRIANGULAR DECOMPOSITION 

Throughout this paper let k be a field of characteristic zero 
and K be the algebraic closure of k. Let k[y] := k[j/i , . . . , y n ] 
be the polynomial ring over the field k in variables yi < 
■ ■ ■ < y n . Let p G k[y] be a non-constant polynomial. The 
greatest variable appearing in p is called the main variable, 
denoted by mvar(p). The integer k such that yu = mvar(p) 
is called the level of p. The separant sep(p) of p w.r.t 
mvar(p), is dp/dmv&r(p). The leading coefficient and the 
leading monomial of p regarded as a univariate polynomial 
in mvar(p) are called the initial and the rank of p; they 
are denoted by init(p) and rank(p) respectively. Let q be 
another polynomial of k[y], we say rank(p) is less than 
rank(q) if mvar(p) < mvar(q), or mdeg(p) < mdeg(g) when 
mvar(p) = mvar(g). 

Let F C k[y] be a finite polynomial set. Denote by (F) 
the ideal it generates in k[y]. Let h be a polynomial in 
k[y], the saturated ideal (F) : h°° of (F) w.r.t h, is the set 
{q G k[y] | 3m G N s.t. h m q G (F)}, which is an ideal in 



k[y]. The polynomial is regular modulo (F) if it is neither 
zero, nor a zerodivisor modulo (F). Denote by V(F) the 
zero set (or algebraic variety) of F in K n . 

Let T C k[y] be a triangular set, that is a set of non-constant 
polynomials with pairwise distinct main variables. We de- 
note by mvar(T) the set of main variables of polynomials in 
T. A variable in y is called algebraic w.r.t. T if it belongs 
to mvar(T), otherwise it is called free w.r.t. T. For a vari- 
able » £ y, we denote by T <v the subsets of T consisting 
of the polynomials t with main variable less than v. Let hr 
be the product of the initials of polynomials in T. We de- 
note by sat(T) the saturated ideal of T: if T is empty then 
sat(T) is defined as the trivial ideal (0), otherwise it is the 
ideal (T) : /iff. The quasi- component W(T) of T is defined 
as V(T) \ V{hr)- Let h G k[y] be a polynomial. Define 
Z(T,h) := W(T)\V(h). 

Let h G k[y] be a polynomial. The iterated resultant of 
h w.r.t. T, denoted by ires(/t,T), is defined as follows: 
(1) if h £ k or all variables in h are free w.r.t. T, then 
ires(/i, T) = h; (2) otherwise, if v is the largest variable of 
h which is algebraic w.r.t. T, then ires(/i, T) = ires(r, T<„) 
where r is the resultant of h and the polynomial in T whose 
main variable is v. Iterated resultants have the following im- 
portant property: the polynomial h is regular modulo sat(T) 
if and only if we have ires(/i, T) 7^ 0. 

We say that the triangular set T is a regular chain if either 
T = or ires(/tT, T) 7^ 0. The pair [T, h] is called a regular 
system if T is a regular chain, and ires(ft, T) 7^ 0. Denote 
by sep(T) the product of all sep(p), where p G T. Then 
T is said to be squarefree if ires(sep(T), T) 7^ 0. A regular 
system rs — [T, h] is said to be squarefree if T is squarefree. 

For a regular system rs = [T,h], the rank of rs, denoted 
by rank(rs), is defined as the set of rank(p) for all p G T. 
Given another regular system rs' = [T' , h'\ with rank(rs) 7^ 
rank(rs'), we say rank(rs) is less than rank(rs') whenever 
the minimal element of the symmetric difference (rank(rs) \ 
rank(rs')) U (rank(rs') \ rank(rs)) belongs to rank(rs). 

A constructive set of K n is any finite union 

{A x \Bx) U ■■■ U (A e \B e ) 

where A\, . . . , A e , B\, . . . , B e are algebraic varieties in K n . 
For any constructible set cs of K n there exist finitely many 
regular systems rsi, . . . , rs m of k[y] such that cs = Z(rsi) U 
••• U Z(rs m ). 

Example 1. Consider the polynomials in k[j/i < yi < y%\ 
Pi = yl +Vi - 1 and p 2 = yxyl - 1. 
We illustrate the previous main notions as follows. 
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The initial yi of p2 is regular modulo (pi). The set T = 
{pi,f>2} is a triangular set. The iterated resultant of y\ and 
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T is 3/1, so T is a regular chain. The pair [T, 3/2] is a regular 
system, since ires(3/ 2 ,T) = y\ — 1. The quasi-component of 
T is the set of points in K 3 such that p 1 = 0, P2 = and 
yi / 0, which is a constructible set. 

We review three important operations Intersect, MakePair- 
wiseDisjoint (MPD) and SymmetricallyMakePairwiseDis- 
joint (SMPD) proposed in [9]. Let rs* = [T*,ft»] be a 
squarefree regular system of k[y] and let p be a polynomial 
of k[y] such that p is regular w.r.t sat(T«). The operation 
lntersect(p, rs») computes a family of squarefree regular sys- 
tems 1Z of k[y] such that 

V(p) n Z(rs*) = UrsenZ{rs), 

and the rank of each rs £ 1Z is less than that of rs*. 

For regular systems [Ti, hi], . . . ,[T e , h e ] in k[y], the function 
MPD returns regular systems [Si, gi], ■ ■ ■ , [Sf, gj] in k[y] s.t. 

Z(Ti,hi) U ■■■ U Z(T e ,h e ) = Z{Sx,g{) U ■■■ U Z{S f ,g f ), 

and for all 1 < i < j < / we have Z(Si,gi) CI Z(Sj,gj) = 0. 

Given a family C = {Ci , . . . , CV} of constructible sets of 
K n , the function SMPD returns a family 7? = {Di, . . . , D s } 
of constructible sets of K" such that D; n Dj = for all 
1 < i < j/' < n, each Dj is a subset of some d, and each d 
can be written as a finite union of some of the -Dj's. Such a 
D is called an intersection-free basis of C. 

3. ZERO SEPARATION 

In this section, we assume n > 2 and we regard the ordered 
variables 3/1 < ■ ■ ■ < y n -i as parameters, denoted by u. Let 
7r u be the projection function which sends a point (u, y~„) 
of K n to the point u of the parameter space K n_1 . Let 
u £ K" -1 . We write 7r~ 1 (u) for the set of all points (u,y n ) 
in K n such that 7r u (u, y n ) = u. 

Let p £ k[u, 3/ n ] be a polynomial of level n, that is, with main 
variable y n . In broad terms, the goal of this section is to de- 
compose the parameter space K n_ into finitely many cells 
such that above each cell the "root structure" of p (num- 
ber of roots, their multiplicity, . . . ) does not change. In 
fact, we make this problem more general by allowing alge- 
braic constraints on the parameter u. After some notations, 
we define in Definition [T] the object to be computed by the 
algorithm devised in this section. It can be seen as a spe- 
cialization of the comprehensive triangular decomposition 
(ctd) to the case where the input system is a regular sys- 
tem and all variables but one are regarded as parameters. 
This algorithm is stated in Section f3. II after two lemmas. 

Notations. Let rs = [T, h] be a regular system of k[u, y n ]. 
If y n does not appear in rs, we denote by Z u (rs) the zero set 
of rs in K n_1 . If y n does not appear in T, we write W a (T) 
for the quasi-component of T in K" _1 . If mvar(/i) = y n 
holds, we denote by coeff(ft) be the set of coefficients of h 
when h is regarded as a polynomial in y n with coefficients 
in k[u] and by Vu(coeff (h)) the variety of coeff (h) in K" _ . 
Finally, if y n is algebraic in T, letting t n be the polynomial 
in T with main variable y n , we write T u = T \ {t n } and 
rs u = [T u ,r], where r = res(h ■ sep(i„),t„) is the resultant 
of h ■ sep(t n ) and t n w.r.t y n . 



Definition 1. Let C be a constructible set of K n_1 . A 
finite set of level n polynomials V C k[u, y n ] separates above 
C if for each a £ C: (I) the initial of any p £ V does not 
vanish at a; (2) the polynomials p(a,y n ) £ K[j/ n ], p £ V, 
are squarefree and coprime. 

Let C be a finite collection of pairwise disjoint constructible 
sets of K n_1 , and, for each C £ C, let Vc C k[u, y n ] be 
a finite set of level n polynomials. Let rs* = [T*,ft»] be a 
regular system of k[u, y n ], where n > 2 and y n is algebraic 
w.r.t T. We say that the family {(C,Vc) C £ C} separates 
Z(rs») if the following conditions hold: 

(1) C is a partition of n xl (Z(rs t )), 

(2) for each C £ C, Vc separates above C, 

(3) Z{rs,) = U CeC U pePc V(p)n ? r- 1 (C). 

More generally, let cs be a constructible set of K n such that 
there exist regular systems rsi, . . . ,rs r of k[u, y n ] whose 
zero sets form a partition of cs and such that y n is algebraic 
w.r.t. the regular chain of rsi, for all 1 < i < r. Then, we 
say that the family {{C,Vc) C £ C} separates cs if C is a 
partition of 7r u (cs) and if for all 1 < i < r there exists a non- 
empty subset d of C and for each C £ d a non-empty subset 
Vc,i C Pc such that {(C,Vc,i) \ C £ d} separates Z(rsi). 
In this case, we have: cs = Ucec U p gp c ^(p) n ""u 

Example 2. Consider the polynomials in k[x > b > a] 

pi = ax 2 — b and pi = ax 2 + 2x + b, 

and the constructible set C = {(a, b) £ K 2 | ab(ab — 1) ^ 
0}. For any point (a, 6) of C, the two polynomials pi(a,6) 
and p2(a,b) of K[k] are squarefree and coprime. So the 
polynomial set {pi,p2} separates above C. 

Consider the regular system rs* = [{pi}, 1] and the con- 
structible sets 

Ci = {(a, b) £ K 2 I ab^ 0} 

C 2 = {(a, 6) £ K 2 I a / & 6 = 0} 

Note that the zero set of rs* is {pi = & o ^ 0}. So the 
family { (Ci,{pi}), (C2,{ax}) } separates Z(rs t ). 

Given two regular systems 

rsi = [{pi}, b] and rs 2 = [{p2,b}, 1]. 
Consider the constructible set 

cs = Z(rsi) U Z(rs 2 ) 

= (V( Pl )\V(ab))u(V(p2,b)\V(a)). 

The family { (Ci,{pi}), (C*2,{p2}) } separates cs. 

Lemma 1. Let p £ k[u,y n ] fee a level n polynomial. Let 
r = res(sep(p) , p) be the resultant of sep(p) and p w.r.t y n . 
Then, the polynomial p(u) o/K[j/ n ] is squarefree and init(p) 
does not vanish at u £ K n_1 , i/ and only if, r(u) 7^ holds. 

Observe that init(p) is a factor of r. So the conclusion follows 
directly from the specialization property of subresultants. 
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Lemma 2. We have the following properties: 

(1) Ify„ does not appear in rs, then ir u (Z(rs)) — Z u (rs). 

(2) If y n does not appear in T and if mvar(h) — y„ holds, 
then we have TT u (Z(rs)) = W U (T) \ V u (coeff (h)). 

(3) If y n is algebraic w.r.tT and if the regular system rs 
is squarefree, then rs u is a squarefree regular system 
o/k[u]; moreover there exists a family 1Z' of squarefree 
regular systems o/k[u,y n ] such that: 

(a) the rank of each rs' £ 1Z' is less than that of rs, 

(b) for each [T',h'] £ 1Z' , y n is algebraic w.r.t T' , 

(b) the zero sets Z(rs'), rs' £ 1Z' and the zero set 
V(t n )nZ(rs u ) are pairwise disjoint, and we have 

(d) Z(rs) = V(t n )nZ(rs u ) U [J rs>e7Z , Z(rs'). 

Proof. Property (1) is clear and proving (2) is routine. 
We prove (3). Since rs is squarefree, using the above nota- 
tions, we have 

ires(r,T) = ires(r,T u ) = ires(/i • sep(t n ),T) / 0. 

This implies that r is regular w.r.t sat(T) and that rs u = 
[Tl^r] is a squarefree regular system of k[u]. Observe now 
that the zero set of rs decomposes in two disjoint parts: 

Z(rs) = (Z(rs) \ V{r)) U {Z{rs) n V(r)) . 

For the first part, we have 

Z(rs)\V(r) = V{t n )nZ(rs u ). 

For the second part, since r is regular w.r.t sat(T), by call- 
ing operation Intersect, we obtain a family 1Z of squarefree 
regular systems of k[u, y n ] such that 

Z(rs)nV{r) = |J Z(rs), 

rj'£R 

where the rank of each rs' £ 1Z is less than that of rs. 
Finally, applying the operation MPD to 1Z we obtain a family 
1Z' satisfying the properties (a), (b), (c) and (d). □ 

3.1 The Algorithm SeparateZe ros 

We present now an algorithm "solving" a regular system in 
the sense of Definition [l] Precise specializations and algo- 
rithm steps follow. 

Calling sequence. SeparateZerosfYs* , u, n) 

Input. A (squarefree) regular system rs* = [T»,/i»] of 
k[u, y n ], where n > 2 and y„ is algebraic w.r.t T„. 

Output. A finite family {(C,Vc) | C £ C}, where C is a 
finite collection of constructible sets of K n_1 , and for each 
C £ C, Vc C k[j/i, . . . , y n ] is a finite set of level n polyno- 
mials, such that {(C, Vc) | C £ C} separates the zero set of 
rs*. (See Definition [T]) 

Step (1). Initialize 1Z = {rs*} and V = 0. 

Step (2). If 1Z — 0, go to Step (3). Otherwise arbitrarily 
choose one regular system rs = [T, h] from 1Z and let 1Z = 



lZ\\rs). Using the above notations, let 72.' be as in Property 
(3) of Lemmai Set V = V U {(rs u , t n )}, set 1Z = 1Z U 1Z' 
and repeat Step (2). 

Comment. Observe that Step (2) will finally terminate since 
each newly added regular system into 1Z has a rank less than 
that of the one removed from 1Z. When Step (2) terminates, 
we obtain a family V of pairs such that 

Z(rs.)= |J V(t n ) Htt' 1 {Z u {rs u )), 

(rs„,l„)eP 

and the union is disjoint. Next, observe that for each pair 
(rs u ,t n ) £ V, the polynomial init(i n ) does not vanish at 
any point of Z u (rs u ), by virtue of Lemma[T] Therefore, the 
union of all Z u (rs u ) is equal to 7r u (^(rs*)). 

Step (3). By means of the operation SMPD we compute an 
intersection- free basis of all Z u (rs u ). Hence we obtain a par- 
tition C of n a (Z(rs*)). Then, for each C £ C we define Vc 
as the set of the polynomials t n such that there exists a reg- 
ular system rs u satisfying (rs u ,t n ) £ V and C C Z u (rs u ). 
Clearly {(C,Vc) \ C £ C} is a valid output. 

Finally, we generalize this algorithm in order to apply it to 
a constructible set represented by regular systems. 

Calling sequence. SeparateZeros({rsi , . . . ,rs r }, u, n) 

Input. Regular systems rsi,...,rs r of k[u, y n ], n > 2, 
whose zero sets are pairwise disjoint and such that y n is 
algebraic w.r.t. the regular chain of rsi, for all 1 < i < r; 
let cs be the constructible set represented by rsi, . . . , rs r . 

Output. A finite family {(C,Vc) \ C £ C}, where C is a 
finite collection of constructible sets of K n_1 , and for each 
C £ C, Vc C k[j/i, . . . ,y n ] is a finite set of level n poly- 
nomials, such that {(C,Vc) | C £ C} separates cs. (See 
Definition [T]) 

Step (1). For each 1 < i < r, call SeparateZeros(rSi, u, n) 
obtaining {(C,Vc) \ C £ Ci} where Ci is a partition of 
ir u (Z(rsi)). 

Step (2). By means of the operation SMPD, compute an 
intersection-free basis T> of the union of the d, for 1 < i < r. 

Step (3). For each D £ T>, let Vd be the union of the Vc 
such that DCC holds. Return {(D, V D ) \ D £ V}. 

4. CYLINDRICAL DECOMPOSITION 

In this section, we propose the notion of an F -invariant 
cylindrical decomposition of K™, generalizing ideas that are 
well-known in the case of real fields. The main algorithm 
and its subroutines for computing such a decomposition are 
stated in three subsections. 

Definition 2. We state the definition by induction on n. 
For n = 1, a cylindrical decomposition of K is a finite collec- 
tion of sets {Dx, . . . , D r +i}, where either r — and Dx = K, 
or r > and there exists r nonconstant coprime squarefree 
polynomials pi, . . . ,p r of k[i/i] such that 

A = {i/ieK|pj(ift) = 0},l<i<r, 
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and D r +i = {j/i £ K | pi(yi) ■ • -p r (yi) 7^ 0}. Note that all 
Di, 1 < i < r + 1 form a partition of K. Now let n > 1, and 
let T>' — {Di, . . . , D s } be any cylindrical decomposition of 
K n_1 . For each Di, let . . . ,f>i,r;}, r« > 0, be a set of 

polynomials which separates above Di. (See Definition [T]) 
If n = 0, set A,i = Di x K. If r-j > 0, set 

A,, = {(a,y n ) £ K" I a e A & p ld (a,y n ) = 0}, 

for 1 < j < n and set 

A,r i+ i = {(a,y n ) 6K"|aefl,& ^nP«,j(a,yn)J # 0}. 

The collection V = {D itj | 1 < i < s, 1 < j < r, + 1} is 
called a cylindrical decomposition of K n . Moreover, we say 
that T> induces T>'. 

Let _F = {/1 , . . . , f 3 } be a finite set of polynomials of k[yi < 
• ■ • < y n ] ■ A cylindrical decomposition T> of K™ is called 
F-invariant if T> is an intersection-free basis of the s + 
1 constructive sets V(fi), 1 < i < s and {y £ K n | 

A(v)---/.(y)#o}. 



Lemma 3. Let rsi, . . . ,rs r +i, with r > 1, be regular sys- 
tems o/k[yi] such that their zero sets form a partition 0/K 1 . 
Then, up to renumbering, there exist polynomials pi, ■ ■ ■ ,p r , 
hi, ... , h r , h r +i £ k[yi] such that rsi = [{pi}, hi] for 1 < i < 
r and rs r +i = [0,h r +i]. Moreover, setting Di = V(pi) for 
1 < i < r and D r +i = {yi £ K | pi(yi) ■ ■ ■ p r {yi) / 0}, the 
sets Di, . . . , -Dr+i form a cylindrical decomposition o/K. 

Proof. Observe that for 1 < i < r we have Z(rSi) = 
V(pi), as hi and pi have no common roots. Since the zero 
sets Z(rsi), . . . , Z{rs r +\) form a partition of K 1 , we must 
have V(h r+ i) — V(pi ■ ■ -p r ). The conclusion follows. □ 

4.1 The Algorithm MakeCyii ndrical 

Calling sequence. MakeCylindrical(7?,, n) 

Input. TZ, a finite family of regular systems such that the 
zero sets Z(rs), for all rs £ 1Z, form a partition of K". 

Output. T>, a cylindrical decomposition of K" such that 
the zero set of each regular system in 1Z is a union of some 
cells in T>. 

Step (1): Base case. If n > 1, go to (2). If 1Z has only 
one element, return T> — K otherwise use the construction 
of Lemma [3] to return a cylindrical decomposition T>. 

Step (2): Initialization. Set to 7Z\,1Zi,lZz the subset of 
7Z consisting of regular systems rs — [T, h] such that, y n is 
algebraic w.r.t T, y n appears in h but not in T, y n does not 
appear in T nor in h, respectively. 

Step (3): Processing TL\. Call SeparateZeros(7?-i, u, n) (see 
Section [3]) obtaining {(C, Vc) I C G Ci} where Ci is a parti- 
tion of 7r u (csi), where csi is the constructible set represented 
by 1Z\. By adding a "1" in each pair, we obtain a collection 
of triples 71 = {(C,Vc, 1) | C € Ci}. 



Step (4): Processing 7?.2- For each rs £ 7?.2, compute the 
projection n u (Z(rs)) by Property (2) of Lemma[2] Set C2 ~ 
{n u (Z{rs)) I rs £ ft 2 } and T 2 = {(C, 0,2) | C £ C 2 }. 

Step (5): Processing 7^.3. For each rs £ 7?-3, compute the 
projection n u (Z(rs)) by Property (1) of Lemma[2] Set C3 = 
{7r u (Z(rs)) I rs £ ^3} and T 3 = {(C, 0, 3) | C £ C 3 }. 

Comment. Since the zero sets of regular systems in 1Z are 
pairwise disjoint, after step (3), (4), (5), we know that the 
element in C3 has no intersection with any element in Ci or 
Ci. Note that it is possible that an element in Ci has inter- 
section with some element of C2. So we need the following 
step to remove the common part between them. 

Step (6): Merging. Set C = C\ UC2 UC3 and T = T 1 UT 2 l>%. 
Note that each element in T is a triple (C,Vc ,%c), with 
C £ C and where 1c is an integer of value 1,2 or 3. By 
means of the operation SMPD, compute an intersection- free 
basis C' of C. For each C £ C' , compute Qc> (resp. Jc 1 ) 
the union of the Vc (resp. Xc) such that C'CC holds. Set 
T' = {(C, Qc,Jc) \CeC}. 

Step (7): Refinement. To each C £ C, apply operation 
MPD to the family of regular systems representing C, so as to 
obtain another family TZc of regular systems representing C 
and whose zero sets are pairwise disjoint. For each rs £ TZc, 
set Vrs = Qc and T rs — Jc- Let 7Z' be the union of the 
TZc, for all C £ C. Set T" = {(Z(rs),V rs ,l rs ) | rs £ TZ'}. 

Comment. Recall that the union of zero sets of the Z(rs), 
for all rs £ TZ equals K". Therefore, it follows from Steps 
(6) and (7), that {Z(rs) \ rs £ TZ'} is a partition of K"~\ 

Step (8): Recursive call. Call MakeCylindrical(7?.', n — 1) to 

compute a cylindrical decomposition T>' of K* 1 such that 
Z(rs), for each rs £ TZ' , is a union of some cells of T>' . For 
each D' £ T>' , observe that there exists a unique rs £ TZ' 
such that D' C Z(rs), so set Pd' = "P rs and T D i — T rs . 
Then, set T" = {(D' ,V D > ,1 D <) \ D' £ V'}. 

Comment. By the comment below Step (5), we know that 
for each triple (D' ,Vo' ,^d') of T'" , the values of J D i can 
only be {1, 2}, {2} or {3}. Next, observe that for each D' £ 
T>' such that X D i = {2} or X D / = {3} holds, we have Vo' — 
0, whereas for each D' £ T>' such that T D i — {1, 2} the set 
Vd 1 is a nonempty finite family of level n polynomials in 
k[j/i, . . . , 2/n] such that Vd' separates above T>' . In Step (9) 
below, we lift the cylindrical decomposition V' of K"~ to 
a cylindrical decomposition T> of K™. 

Step (9): Lifting. Initialize T> to the empty set. For each 
D' £ V such that 1 D , = {2} or I D , = {3} holds, let V := 
V U {D' x K}. For each D' £ V such that = {1, 2} 
holds, let I? = V U {D P }, where 

D p = {(q, y„) £ K" \ a E D' and p(a, y n ) = 0}, 

for each p £ Vo 1 an d let I? = T) U {-D*}, where 

D* = {(a,y„)£K"|a£D'& j p(et,y n ) J / 0}, 

Finally, return X>. The correctness of the algorithm follows 
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from all the comments and Definition [2] 

4.2 The Algorithm InitialPartition 
Calling sequence. I n itial Partition (_P, n) 

Input. F — {/i, . . . , / s }, a finite subset of k[yi < • • • < y n ]. 

Output. A family TZ of regular systems, the zero sets of 
which form an intersection-free basis of the s+1 constructible 
sets n/i)> • • • . V(f a ) and {y G | (]JU Mv)) + °>- 

Step (1): Let B = SMPDfV (,fi), . . . , V{f 3 )) be an intersec- 
tion free basis of the s constructible sets V(fi), . . . , V(f s ). 
For each element B of B, we apply operation MPD to the 
family of regular systems representing B to compute another 
family TZb of squarefree regular systems such that the zero 
sets of regular systems in TZb are pairwise disjoint and their 
union is B. Let 7Z be the union of all TZb, B G B. Clearly 
the set {Z(rs) \ rs G 7Z} is an intersection-free basis of the 
s constructible sets V(fi), . . . , V(f s ). 

Step (2): Let / = H f . eF fi and rs, = [0,/]. Set TZ = 
1Z U {rs*}. Obviously TZ is the valid output. 

4.3 The Algorithm Cyli nd rica I Decom pose 

Calling sequence. CylindricalDecompose(7 ? , n) 

Input. F, a finite subset of k[j/i < ■ • • < y n ]. 

Output, an F-invariant cylindrical decomposition of K n . 

Step (1): If n > 1, go to step (2). Otherwise let {pi, . . . ,p r }, 
r > 0, be the set of irreducible divisors of non-constant 
elements of F. If r = 0, set T> — K and exit. Otherwise set 

A = {yi e K | jJiGrn) = 0}, l < i < r, 

and D r +i = {j/i 6 K | pi(yi) ■ ■ ■ p r (yi) 0}. Clearly 
T> — {Di |l<i<r + l}isan F- invariant cylindrical 
decomposition of K. 

Step (2): Let TZ be the output of lnitialPartition(_F, n). 

Step (3): Call algorithm MakeCylindrical(7?,, n), to compute 
a cylindrical decomposition T> of K n such that the zero 
set of each regular system in TZ is a union of some cells 
in T>. Clearly, T> is an intersection-free basis of the set 
{Z(rs) | rs G TZ}, which implies T> is an intersection-free 
basis of the s + 1 constructible sets V(fi), . . . ,V(f a ) and 
{y € K" | (TK=i /*(»)) * °>- Therefore, © is an F- 
invariant cylindrical decomposition of K n . 

5. CYLINDRICAL ALGEBRAIC DECOM- 
POSITION 

In this section, we show how to compute a CAD of R n from 
a cylindrical decomposition of C n . This section starts with 
reviewing basic notions for CAD pQ. A theorem (Theorem[T} 
due to Collins [TT] is then reviewed, where the relation be- 
tween complex and real roots of a polynomial with real co- 
efficients is shown. The bridge from cylindrical decomposi- 
tion to CAD is built in Corollary [T] which can be directly 
obtained from Collins' theorem. The main algorithm CAD 
and its subroutines are stated in four subsections. 



A semi- algebraic set !4 4 of R n is a subset of R n which can 
be written as a finite union of sets of the form: 

{y G R n | V/ G F, f(y) = and V 5 G G,g(y) > 0}, 

where both F and G are finite subsets of the polynomial 
ring R[y u . . . ,y n ]. 

Given an n-dimensional real space R n , a nonempty con- 
nected subset of K" is called a region. For any subset S 
of R™, a decomposition of S is a finite collection of disjoint 
regions whose union is S. For a region R, the cylinder over 
R, written Z(R), is R x R 1 . Let /i <■■■</,-, r > be con- 
tinuous, real- valued functions defined on R. Let fo = — oo 
and fr+i ~ +oo. For any fi, 1 < i < r, we call the set of 
points {(a, fi(a)) | a G R} the fi-section of Z(R). For any 
two functions fi, fi+i, < i < r, the set of points (a, 6), 
where a ranges over R and fi(a) < b < fi+i(a), is called 
the (fi, fi+i) -sector of Z(R). All the sections and sectors of 
Z(R) can be ordered as 

(f0, fl) < /l < • • • < /r < (fr, fr+l). 

Clearly they form a decomposition of Z(R), which is called 
a stack over R. 

A decomposition £ of R" is cylindrical if either (1) n = 1 
and £ is a stack over R°, or (2) n > 1, and there is a cylin- 
drical decomposition £' of R n_1 such that for each region 
R in £' , some subset of £ is a stack over R. Moreover, We 
say that £ induces £' . A decomposition is algebraic if each 
of its regions is a semi-algebraic set. A cylindrical algebraic 
decomposition of R n is a decomposition which is both cylin- 
drical and algebraic. 

Let p be a polynomial of R[j/i, . . . , y n ], and let S be a subset 
of R n . The polynomial p is invariant on S (and S is p- 
invariant) , if the sign of p(a) does not change when a ranges 
over S. Let F G R[j/i, . . . ,y n ] be a finite polynomial set. 
We say S is F-invariant if each p G F is invariant on S. A 
cylindrical algebraic decomposition £ is F-invariant if F is 
invariant on each region of £ . 

Let p be a polynomial of R[j/i, . . . , y n ], and let R be a region 
in R n_1 . p is delineable on 7? if the real zeros of p define 
continuous real- valued functions 9i, . . . , 9 S such that, for all 
a G R, 9i(a) < ■■■ < 9 s (a), and for each there is an 
integer m, such that mi is the multiplicity of the root di(a) 
of p(a,y n ). Note that if k = 0, V(p) has no intersection 
with . Clearly when p is delineable on R, its real zeros 
naturally determine a stack over R. 

Let £ be a CAD of R n . As suggested in [1], each region e G £ 
can be represented by a pair (7, S), where I is the index of 
e and S is a sample point for e. The index 7 and the sample 
point S of e are defined as follows. If n = 1, let 

ci < e 2 < ■ ■ ■ < e 2m < c 2m +i, m > 

be the elements of £ . For each ej, the index of ej is defined 
as (i). For each e;, its sample point is any algebraic point 
belonging to ei. Let £' be the CAD of R" induced by £ . 
Suppose that region indices and sample points have been 
defined for £ ' . Let 

e;,i < e<,2 < • • • < ei,2 mj < ej^mi+i) > 
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be the elements of £ which form a stack over the region e< of 
£' . Let (ii, . . . ,i„-i) be the index of ej. Then the index of 
eij is defined as (ii,.. . ,i n -i,j). Let 5" be a sample point 
of e». Then the sample point of e;,j is an algebraic point 
belonging to dj such that its first n — 1 coordinates are the 
same as that of S' . 



Theorem 1 (Collins). £ei p fee a polynomial of ring 
< ••• < y n ] R be a region o/R 71-1 . If init(p) ^ 
on R and the number of distinct complex roots of p is 
invariant on R, then p is delineable on R. 

Corollary 1. Let F = {pi, ■ ■ ■ ,Pr} be a finite set of 
polynomials in R[j/i < ■ ■ ■ < y n ] of level n. Let R be a 
region o/R n_1 . Assume that for every a £ R, (1) the initial 
of each pi does not vanish at a; (2) all Pi(a, y n ), 1 < i < r, 
as polynomials ofM.[y n ], are squarefree and coprime. Then 
each pi is delineable on R and the sections of Z(R) belonging 
to different pi and pj are disjoint. 

Let R and F be defined as in the above corollary. Then 
clearly the real roots of all p £ F are continuous functions 
on R and they together determine a stack over ii. The 
algorithm GenerateStack, described in Section f5. 2 1 is a direct 
application of the above corollary. 

5.1 Real Root Isolation 

Let a = (fin, ■ ■ ■ ,a n ) be an algebraic point of R n . Each cti 
as an algebraic number is a zero of a nonconstant squarefree 
polynomial t;(i/i) of Q[j/i]. Let T be the set of all ti(j/j). 
Clearly T is a zero dimensional squarefree regular chain of 
Q[y]. On the other hand, if T is a zero-dimensional regular 
chain of Q[y], any real zero of T is an algebraic point of R n . 
Therefore any algebraic point a of R n can be represented by 
a pair (T, L) , where T is a zero-dimensional squarefree reg- 
ular chain of Q[y] such that T(a) = and L is an isolating 
cube containing a but not other zeros of T. The pair (T, L) 
is called a regular chain representation of a, which will be 
used to represent a sample point of CAD. 

Next we provide the specification of an algorithm called Iso- 
lateZeros for isolating real zeros of univariate polynomials 
with real algebraic number coefficients. It is a subroutine of 
the algorithm NREALZERO proposed in [34] for isolating the 
real roots of a zero-dimensional regular chain. 

Calling sequence. lsolateZeros(a' n ~ 1 ' , F, n) 

Input, a'™" 1 ' is a point of R n_1 , n > 1, with a regular 
chain representation (T',L'). If n = 1, T' = and L' = 
0. F — {pi, . . . ,p r } is a list of non-constant polynomials 
of Q[j/i,-- - ,y n ] of level n satisfying that (1) for pi 6 F, 
T' U {pi} is a squarefree regular chain of Q[yi, . . . , y n \\ (2) 
all piia^-V^n), 1 < i < r, as polynomials of R[y n ], are 
squarefree and coprime. 

Output. A pair (TV, u). Let p = Ui=i Pi- N = (Ni, . . . , N m ) 
is a list of intervals with rational endpoints with Ni < ■ • ■ < 
N m such that each Nj contains exactly one real zero of 
p(a ( - n ~ 1 \y n ). v = (yi, . . . , u m ) is list of integers, where 



1 < fi < r, such that the zero of p(a' n x \yn) in Nj is a 
zero of (a (n_1) , y n ). 

5.2 The Algorithm GenerateStack 
Calling sequence. GenerateStack(e', F, n) 

Input, e is a region of a CAD £' of R n_1 , n > 1, and 
e' is represented by its index /' and its sample point S' . 
Let (T',L') be the regular chain representation of S". If 
n = 1, I',T',L' = 0. F is a finite set of polynomials in 
Q[yi, ■ ■ ■ ,y n ] of level n. The region e' and the polynomial 
set F satisfy the conditions specified in Corollary [T] 

Output. A stack S over e'. 

Step (1). If F = 0, go to step (2). Otherwise call algorithm 
lsolateZeros(<5", F, n) to isolate the real roots of polynomials 
in F w.r.t y n at the sample point S' of e'. Let (N, v) be the 
output. If TV 7^ 0, go to step (3). 

Step (2). Let I = (/', 1). Let T = T' U {y n }, L = L 1 x [0, 0], 
S = (T, L) and return S = ((/, 5")). 

Step (3). Let Ni = toi,6i],...,iV m = [a m ,b m ], m > be 
the elements of N. For 1 < i < 2m + 1, set Ii = (I',i). 
Let Si be the greatest integer less than oi. Let S2 m +i be 
the smallest integer greater than b m . For 1 < i < m — 1, 
let s 2i+ i = For < i < m, Let T 24+ i = T' U 

{y n - «2i+i}, L 2 i+i — L' x [s2i+i, sai+i] and set Szi+i = 
(T 2 i + i,L 2 i+i). For 1 < i < m, let T 2i = T 1 Up„„ L 2i = 
L' x Nt and set &, = (T 2 i, L 2i ). Finally, set S be the list of 
all (L, Si), 1 < i < 2m + 1. Then 5 is the stack over e'. 

5.3 The Algorithm MakeSemiAlgebraic 
Calling sequence. MakeSemiAlgebraic(I?, n) 

Input. I? is a cylindrical decomposition of C n , n > 1. 

Output. A CAD f of R n such that, for each element D of 
T>, the set D n R n is a union of some regions in £. 

Step (1). If n > 1 goto (2). Otherwise let Di, ... ,D r ,D r+1 , 
r > be the elements of D. For each 1 < i < r, let Pi be 
the polynomial such that Di = {yi | pi(yi) = 0}. Let £ be 
the output of GenerateStack(0, {pi, . . . ,p r }, 1). Clearly £ is 
a CAD of R 1 . 

Step (2). Let 23' be the cylindrical decomposition of C' 1 " 1 
induced by T>. Call MakeSemiAlgebraic recursively to com- 
pute a CAD £' of R n_1 . 

Step (3). In this step we lift the CAD £' of R n_1 to £. 
Initialize £ = ( ). For each region e of £' , let D be the 
cell of V such that e' C D' fl I". Let D\, . . . , D r ,D r+ i, 
r > be the cells of V such that D' x C = Uj±jDj. 
For each 1 < j • < r, let be the polynomial such that 
Dj = {(a, y n ) \ a £ D' &c Pj(a,y n ) = 0}. Add output of 
GenerateStack(e', {pi, . . . ,p r },n) into £ . Clearly £ is a CAD 
of R™ and for each D £T>, the set D n R n is a union of some 
regions in £. 

5.4 The Algorithm cad 
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Calling sequence. CAD(_F, n) 

Input. F is a finite subset of Q[yi < ■ ■ ■ < y n ], n > 1. 
Output. An F-invariant CAD £ of R n . 

Step (1). Let T> — CylindricalDecompose(F, n) be an F- 
invariant cylindrical decomposition of C n . 

Step (2). Call algorithm MakeSemiAlgebraic to compute a 
CAD £ of R n such that, for each element D of V, the set Dn 
R™ is a union of some regions in £ . Since T> is an intersection- 
free basis of the s + 1 constructible sets Vc(/i), • • ■ , Vc(fs) 
and {y £ C" | (Tii=i fi(y)) / 0}, £ is an intersection-free 
basis of the s + 1 semi-algebraic sets ■ • ■ , Vr(/ s ) and 

(j £ 1" (n<=iA(f)) ^ °i- Note that each ele ment in 
£ is connected. Therefore £ is an F-invariant cylindrical 
algebraic decomposition of R™ . 

6. EXAMPLES AND EXPERIMENTATION 
6.1 An Example 

Let us illustrate our method by a simple and classical exam- 
ple. Consider the parametric parabola p = ax 2 + bx + c. Set 
the order of variables asx>c>b>a. The first step Ini- 
tialPartition generates four regular systems, whose zero sets 
form a partition of C 4 . 




r 3 := 



ax + bx + c — 



r 2 := 




:= { ax 2 + bx + c ^ 



Next we trace the algorithm MakeCylindrical. Initialize the 
sets 72.1 := {r2,rs}, 1Z 2 := {^4} and 72.3 := {fi}- Since x 
appears in the equations of r 2 and r 3 , SeparateZeros(7?.i) is 
called to obtain a family of pairs 

{(C 1 ,{t}),(C 2 ,{p}),(C 3 ,{'l})}, 
defined as follows, which separates Z(r 2 ) U Z(rs). 



Ci 
C 2 

c. 



{a = 0,6/0} -> {t} 

{a(4ac-b 2 )^0} -» {p} 
{4ac - b 2 = 0, a ± 0} -> {q} 



{fcr. + c} 

{aa; 2 + bx + c} 

{2ax + 6} 



The projection of Z(r±) is the values such that a, b, c do 
not vanish simultaneously, denoted by C4. The projection 
of Z(ri) is the set {a — b = c = 0}, denoted by C5. 

Note that Ci,C 2 , C3 are all subsets of C4. In the Merging 
step, by calling SMPD, we get another set Ce := {a = b = 
0, c / 0} such that Ci,C2,C3,C5 and C6 are pairwise dis- 
joint and their union is C 3 . Moreover, for each d, there is 
a family of polynomials and indices associated to it. 



Ci 


c 2 


c 3 


c 5 


C*6 


{t} 


M 


{1} 








{1,2} 


{1,2} 


{1,2} 


{3} 


{2} 



Since each d is already the zero set of some regular system, 
MakeCylindrical({Ci, C 2 , C 3 , C* 5 , C* 6 }, 3) 



is called recursively to compute a cylindrical decomposition 
of C 3 . By the Lifting step, we finally obtain a p- invariant 
cylindrical decomposition of C 4 . Let r = 4ac — b 2 , the de- 
composition can be described by the following tree, 
root 



a =j£ 

\ 

C 



a = 

/\ 

b = 67^0 

\ \ A /\ 

To compute a p-invariant CAD of R 4 from the above tree 
is straightforward. Starting from the root, one first obtains 
the trivial zero of a, which decomposes a / into two 
connected cells a < and a > 0. The real line is thus di- 
vided into three parts. For each part, one then substitutes 
its sample point into its children which are equations, from 
where one can determine the sample points for the children 
which are inequations. Continuing in this manner, one fi- 
nally obtains a CAD of R 4 with 27 cells. The number of cells 
is the same as that obtained in [6]. In fact, it is the minimal 
number of cells one can obtain for a p- invariant CAD of R 4 . 

6.2 Experimental Results 

In this section, we present experimental results obtained 
with an implementation of the algorithms presented in this 
paper. Our code is in Maple 12 running on a computer with 
Intel Core 2 Quad CPU (2.40GHz) and 3.0GB total memory. 
The test examples, listed in appendix for the reader's con- 
venience, are taken from diverse papers [161 [T] 1131 1251 El 1141 
112] on CAD. The time-out for a test run is set to 2 hours. 

In Table 1, we show the total computation time of CAD 
and the time spent on three main phases of it, which are 
InitialPartition (Partition for short), MakeCylindrical (M.C. for 
short) and MakeSemiAlgebraic (M.S. A. for short). We also 
report the number of elements (A/ H ) in the CAD. Aborted 
computations due to time-out are marked with From 
the table, one can see that, except examples 14 and 16, the 
steps of the algorithm dedicated to computations in complex 
space dominate the step taking place in the real space. 

In Table 2, we show the total computation time of the al- 
gorithm Cylindrical Decompose (CD. for short) and the time 
spent on three main operations of it, which are respectively 
SeparateZeros(Separate for short), MPD and SMPD. We can 
see that the cost of algorithm Cylindrical Decompose is domi- 
nated by SMPD. The number of elements (Nc) in the cylin- 
drical decomposition of C n is also reported. 

The data reported in two tables shows that SMPD is the 
dominant operation, which computes intensively CCDs of 
polynomials modulo regular chains. This suggests that the 
modular methods and efficient implementation techniques 
in [15] [24] [22] (use of FFT-based polynomial arithmetic, . . . ) 
have a large potential for improving the implementation of 
our CAD algorithm. 
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Sys Partition M.C. M.S.A. Total N R 



1 


0.024 


0.096 


0.024 


0.144 


27 


2 


1.184 


2.856 


1.048 


5.088 


895 


3 


0.004 


7.512 


0.704 


8.220 


233 


4 


0.264 


1.368 


1.080 


2.716 


421 


5 


0.016 


0.052 


0.116 


0.184 


55 


6 


0.108 


0.156 


0.120 


0.384 


41 


7 


2.704 


3.600 


1.360 


7.664 


893 


8 


0.380 


1.608 


1.196 


3.184 


365 


9 


0.288 


0.532 


0.264 


1.084 


209 


10 


5.668 


48.079 


18.833 


72.640 


3677 


11 


0.252 


1.192 


0.620 


2.068 


563 


12 


2.664 


135.028 


88.142 


225.862 


20143 


13 


10.576 


35.846 


6.905 


53.335 


4949 


14 


5.728 


71.760 


2520.354 


2597.878 


27547 


15 


690.731 


2513.817 


299.250 


3503.954 


66675 


16 


895.435 


2064.469 








17 


0.052 










18 












Table 1 Timing (s) and number of cells for CAD 


Sys 


Separate 


MPD 


SMPD 


Total 


N C 


1 


0.020 


0.012 


0.084 


0.156 


8 


2 


0.508 


0.252 


2.268 


4.052 


63 


3 


3.856 


0.836 


2.460 


7.880 


24 


4 


0.280 


0.088 


1.036 


1.648 


65 


5 


0.032 


0.008 


0.012 


0.064 


7 


6 


0.036 


0.012 


0.092 


0.268 


13 


7 


1.100 


0.652 


2.416 


6.320 


58 


8 


0.536 


0.144 


1.040 


2.008 


55 


9 


0.120 


0.032 


0.384 


0.816 


26 


10 


3.204 


0.756 


49.031 


54.119 


594 


11 


0.128 


0.032 


0.960 


1.416 


49 


12 


8.508 


2.024 


125.104 


138.188 


856 


13 


2.040 


1.784 


42.578 


47.002 


407 


14 


5.741 


2.092 


64.875 


76.956 


983 


15 


83.469 


62.736 


3066.071 


3232.073 


2974 


16 


66.516 


377.664 


2501.947 


2959.904 


5877 



Table 2 Timing (s) and number of cells for CD. 



7. CONCLUSION 

We have presented a new approach for computing cylindrical 
algebraic decompositions. Our main motivation is to under- 
stand the relations between CADs and triangular decompo- 
sitions, studying how the efficient techniques developed for 
the latter ones can benefit to the former ones. 

Our method can be applied for solving QE problems directly. 
However, to solve practical problems efficiently, our method 
needs to be equipped with existing techniques, like partially 
built CADs, for utilizing the specific feature of input prob- 
lems. Such issues will be addressed in a future paper. 
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12. X-axis ellipse problem 

r.2/ \2 i 2 2 2i2 

{b {x — c) + a y — a b , 

x 2 + y 2 — l},y > x > b > c > a. 



13. Davenport and Heintz 
{a — d , b — c , a — c, b — 1, a 2 



b},a > b > c > . 



14. Hong-90 

{r + s + t, rs + st + tr - 

t>s>r>b>a. 



a. rst — b}. 



15. SolotarefT-3 

{r, r — 1, ii + 1, it — v, v — 1, 

3u 2 + 2ru — a, 3v 2 + 2rv — a, 

u ?> + ru 2 — au + a — r — 1, 

v 3 + rv 2 — av — 2b — a + r + 1}, 

b>u>v>r~>a. 

16. Collision problem 
{tJ*-6, %t- 10,x-i|t+ 1, 

x- igt- l,s/- iZt + 9,1/- %t + 7, 
(x - t) 2 + y 2 - l},t > x > y. 

17. McCallum trivariate random polynomial 

{(y - l)z 4 + xz 3 + x(l - y)z 2 + (y - x - \)z + y}, 
z > y > x. 



18. Ellipse problem 
{6 2 (x-c) 2 + o 2 (y-d) 
y>x>d>c>b>a. 



■ a 2 b 2 ,a, b,x 2 + y 2 - 1} , 



APPENDIX 

1. Parametric parabola 

{ax 2 + bx + c}, x > c > b > a. 

2. Whitney umbrella 

{x — uv, y — v, z — u 2 }, v>u>z>y>x. 

3. Quartic 

{x 4 + px 2 + qx + r}, x > p > q > r. 

4. Sphere and catastrophe 

{z 2 + y 2 + x 2 - 1, z 3 + xz + y}, x > y > z. 

5. Arnon-84 

{y 4 - 2y 3 +y 2 - 3x 2 y + 2a: 4 } , y > x. 

6. Arnon-84-2 

{144y 2 + 96a; 2 y + 9a; 4 + 105a: 2 + 70a; - 98, 
xy 2 + 6a:y + x ?> + 9a;}, y > x. 

7. A real implicitization problem 

{x — uv.y — uv 2 , z — u 2 },v > u > z > y > x. 

8. Ball and circular cylinder 

{a: 2 + y 2 + z 2 - l,x 2 + (y + z - 2) 2 - 1}, z > y > x. 

9. Termination of term rewrite system 

{x -r,y - r,x 2 {l + 2y) 2 - y 2 {l + 2a; 2 )}, r > x > y. 

10. Collins and Johnson 

{3a 2 r + 3b 2 - 2ar - a 2 - b 2 , 

3a 2 r + 3b 2 r - 4ar + r - 2a 2 - 26 2 + 2a, 

a - 1/2, b, r, r - 1}, r > a > b. 

11. Range of lower bounds 

{a, az 2 + bz + c, aa; 2 + fra; + c — y}, 
z>c>fc>a>a;>y. 
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